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Abstract 

In this paper we consider the minimum time population transfer problem for a two level quantum 
system driven by two external fields with bounded amplitude. The controls are modeled as real 
functions and we do not use the Rotating Wave Approximation. After projection on the Bloch 
sphere, we tackle the time-optimal control problem with techniques of optimal synthesis on 2- 
D manifolds. Based on the Pontryagin Maximum Principle, we characterize a restricted set of 
candidate optimal trajectories. Properties on this set, crucial for complete optimal synthesis, are 
illustrated by numerical simulations. Furthermore, when the two controls have the same bound 
and this bound is small with respect to the difference of the two energy levels, we get a complete 
optimal synthesis up to a small neighborhood of the antipodal point of the starting point. 
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I. INTRODUCTION 



In this paper we apply techniques of optimal synthesis on 2-D manifolds to the population 
transfer problem for a two- level quantum system (e.g. a spin 1/2 particle) driven by two 
external fields. Two-level systems are the simplest quantum mechanical models interesting 
for applications, see for instance (2J [H)] . The dynamics is governed by the time dependent 
Schrodinger equation (in a system of units such that h — 1): 



where E is a real number (±E represent the two energy levels of the system). The controls 
(^i(-), f^O)) assumed to be real valued and different from zero only in a fixed interval, 
represent external pulsed fields. The Hamiltonian without external fields, i.e., the matrix 
diag(— E, E), is called drift term. 

The goal is to steer the system from the first level (i.e. \ipi\ 2 = 1) to any other target 
state in minimal-time and with controls of bounded amplitude, 



where T is the transfer time, Mi and M 2 are two positive real constants representing maxi- 
mum available amplitudes for the control fields. The most interesting target state is of course 
the second level (i.e. \ip2\ 2 = 1)- 

Remark 1. The two real controls represent two independent fields acting in two orthogonal 
directions. They do not come from the use of the Rotating Wave Approximation close to 
the Bohr frequency of the system as it often happens in problems with two controls. Each 
field acts independently and has its own bound on the amplitude. As a consequence, the use 
of the interaction picture does not permit to eliminate the drift term. More precisely, the 
system would be driftless in the interaction picture, but with a control set depending explicitly 
on time (and not anymore of the form < Mi with Mi constant,). 




(1) 



where = (^i(-), M')) T ■ [0,^] >-> C 2 satisfies £ 2 =1 |^j(i)| 2 = 1, and 




Qi(t)\ <M h i = l,2 for every t e [0,T], 
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The time optimal problem for two level quantum system with one bounded real control 
was studied in [7J. For the same problem with unbounded control, see [12J. The minimum 
energy problem with one unbounded control was addressed in [TTJ. For the minimum energy 
problem with two unbounded controls see [6j [TT] . Regarding optimal control problems for 
two-level dissipative systems, see [31 H]. Surprisingly the time optimal problem with two 
bounded real controls for closed two-level systems has not yet been studied. This is a 
relevant problem in NMR, see [U [13] and references therein. 

It is standard to eliminate global phase by projecting the system on a two dimensional 
real sphere S 2 (called the Bloch Sphere) by means of a Hopf map [7j. After setting Ui(t) = 
tti(t)/Mi, the controlled Schrodinger equation ([T]) becomes a two-input control-affine system 
on the sphere S 2 : 

x = Fx + U\G\X + u 2 G 2 x, \ui(t)\ < 1, 



(3) 



where x := (xi, x 2 , ^3) T £ 
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with a := arctan( v /M 1 2 + M$/E), (5 := arctan(M!/M 2 ), and k := 2^E 2 + Mf + M$. 

Normalizations. To simplify the notation, we normalize k = 1. This normalization 

corresponds to a time re-parameterization. More precisely, if T is the minimum time to steer 

the state one to a target state for the system with k = 1, the corresponding minimum time 

T 

for the original system is . 

2v/£ 2 + M 2 + M 2 

Assumptions. Two types of assumptions on the parameters a and are used in this 
paper: 



(Al) < a < tt/4 and < /3 < tt/4. 
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(A2) a small and /3 = 7r/4. 



Assumption (Al) is used in Sec. III. Assumption (A2) is used in Sec. IV and Sec. [VJ Note 



that in (Al) it is not restrictive to assume < < 7r/4, as the controls U\ and 112 play a 
symmetric role; in (A2) (3 = tt/A corresponds to the case where the two controls have the 
same bound (Mi = M 2 ). 

Remark 2. Roughly speaking, the parameter a measures the relative strength of the control 
fields compared to the static one. The parameter (5 characterizes the relative strength between 
the two control fields. In spin experiments, the static field represented by the drift term is 
many orders of magnitude larger than the radio-frequency control fields JT^j , Chap. 10]. 
Therefore, the most relevant case corresponds to small a. 

The vector fields Fx, G\X, and G 2 x describe rotations respectively around the axes x 3 , Xi, 
and £2- The state one which corresponds to the lowest energy level is represented by the point 
N := (0, 0, 1) (called north pole) and the state two which corresponds to the highest energy 
level is represented by the point S := (0,0,-1) (called south pole). The optimal control 
problem we are interested in is to connect the north pole to any other fixed state in minimum 
time. The most important final state is of course the south pole. In the case of a spin 1/2 
particle the later case corresponds to a complete spin flip. As usual we assume control ««(•) 
to be a measurable function satisfying ]««(£) | < 1 almost everywhere. The corresponding 
trajectory is a Lipschitz continuous function x(-) satisfying ^ almost everywhere. Since ^ 
is controllable, and the set of velocities V(x) := {Fx + U\G\X + U2G1X, \u\\ < 1, \ui\ < 1} 
is compact and convex, solutions to the time-optimal control problem exist. [H Chap. 10]. 
By solution we mean an optimal synthesis, i.e., the collection of time-optimal trajectories 
starting from the north pole: 

{7x| 7x is time optimal between N and x}x<=s 2 - 

Optimal syntheses are considered as the right concept of solutions for optimal control prob- 
lems, see [IB] . One of the most important tools for the construction of optimal synthesis is 
the Pontryagin Maximum Principle (PMP for short, see [T7], [TJ Chap. 12]). It is a first 



order necessary condition for optimality that allows us to restrict the set of candidate op- 
timal trajectories. One then needs to select the optimal ones from this set. In general the 
selection step is the most difficult [HI EES] • 

For arbitrary values of a and /3 satisfying Assumption (Al), we are mainly concerned with 
the fist step, i.e., the construction of a restricted set of candidate time-optimal trajectories. 
The most difficult task is to analyze the role of singular trajectories. For small a and /3 = 7r/4 
(i.e with controls bounded on the square: \Qi\ < M, i = 1,2), we complete the time optimal 
synthesis up to a neighborhood of order a of the south pole. More precisely, the optimal 
synthesis is composed of four families of trajectories starting from N with (u± — 1, 112 — 1), 
or {ui — l,u 2 — — I)? or ( u i — — 1) M 2 = — 1), or {ui = — l,w 2 = 1), an d switching for the 
first time at s with s G [0, s max ] and then every v (s) until reaching a neighborhood of S. 
Here switching means that one of the two controls switches from +1 to —1 or vice versa. 
The expressions of s max and v(s) are given by 



arccos 



sin 2 a 



v(s) — arccos 



1 + cos 2 a / 
d-A(s)-B(s)-C(s) 



e- A(s) + B(s) 

where A(s) := 8 cos a sin 2 a sin(s), B(s) = 2 sin 2 2a cos(s), C(s) := 4 sin 4 a cos(2s), d 



sin 2 2a, and e := 5 + 2 cos 2a + cos 4a. See Proposition HI. 2 and Corollary |III.3 for more 
detail. An image of the time optimal synthesis is presented in Fig. [TJ where the non- 
intersecting four-snake structure of the four families of optimal trajectories is preserved 
outside a neighborhood of the south pole. The colored curves called switching curves are 
the locations where a switching occurs. Note that the switching curves are located close to 
the two great circles passing through the north pole and containing X\— or X2— axis; the 
endpoints of each switching curve are located exactly on these great circles. We deduce 
from the optimal synthesis that for a given target state, the optimal trajectories are bang- 
bang, and the corresponding optimal controls are periodic on all interior bang arcs. In other 
words, u\ and 112 are periodic except on the first and last pieces. The first and the last 
switching times need to be computed numerically depending on the target. These optimal 
trajectories are more complicated than the ones with the two controls bounded on the circle 



X1 



FIG. 1. Optimal synthesis for a = 0.25 and j3 = ir/4. Black curves are optimal trajectories and 
colored curves are switching curves. 

(i.e. with \JVl\ + < M, see [61 E]), but permit faster transfer times. We show in Fig. [2] 
a time-optimal control that steers ^ from iV to 5* with a = 0.25. 

Based on the optimal synthesis up to a neighborhood of S, we also propose two families 
of simple suboptimal controls which both allow a transfer from N to S faster than the 
optimal controls bounded on the circle. More specifically, consider the sequence of controls 
(u\ = 1,1*2 = -1) — > ( u i = — 1 ? ^2 = —1) — > ( u i = —l, u 2 — 1) — > (ui — l,u 2 — 1), or any of 
the three other cyclic permutations of it, where each pair of controls lasts for a duration equal 
to 7r/2. We show that successively applying this sequence steers system ^ from iV to a point 
close to S with an error of order a. Furthermore, a similar but slightly non-saturate sequence 
(ui = 7,«2 = -7) -> (ui = -7,112 = -7) -> (wi = — 7,«2 = 7) -> (wi = 7,m 2 = 7), where 
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FIG. 2. For a = 0.25 and f3 = tt/4, optimal control (u\,U2) (left) and corresponding optimal 
trajectory reaching the south pole (right). 



7 is an explicitly computable positive constant slightly smaller than 1, will bring ([3]) from 
N to S exactly. Both strategies realize a transfer time close to the optimal one without any 
computation. 



The paper is organized as follows. In Section [IT], we derive basic facts of optimal syntheses 
on 2-D manifolds for control-affine systems with two bounded controls. The section is self- 
contained, and has its own value beyond the optimal control problem considered in this 
paper. Based on these results, we present in Section [HI] a restricted set of candidate optimal 



trajectories for the case with < a < ir/4 and < f3 < it/ 4. In Section IV, we complete the 
time optimal synthesis up to a neighborhood of order a of the south pole for small a and 
f3 = 7r/4. Further, we derive in Section [V] two simple suboptimal strategies, and compare 
them with the optimal strategy for controls bounded on the circle. Finally, we gather in 
Appendices [A] [B] and [C] technical proofs and computational lemmas. 



II. OPTIMAL SYNTHESES ON 2-D MANIFOLDS WITH TWO BOUNDED CON- 
TROLS 



In this section, we introduce important definitions and develop basic facts about optimal 
syntheses on 2-D manifolds for control-affme systems with two bounded controls. We use 
ideas similar to those used by Sussmann, Bressan, Piccoli and the first author in [5j [8j El HEJ 
[19]. This section is written to be as self-contained as possible. 

A. Basic Definitions and PMP 

We focus on the following problem: 
(P) Consider the control system 

x = F(x) +u 1 Gi(x) + u 2 G 2 (x), (4) 

where x G M, \ui\ < 1, i = 1,2. We make the following assumption: 

(HO): M is a smooth 2-D manifold. The vector fields F, G\ and G 2 are C°° , and the control 
system is complete on M . 

The goal is to reach every point of M in minimum time from a source M in which is assumed 
to be a smooth submanifold of M , possibly with a smooth boundary. 

In the following we use the notation u := (ui,u 2 ), and x := (xi,x 2 ) in a local chart. 

Definition 1. A control for the system is a measurable function u(-) = (ui(-) , u 2 (-)) ■ 
[di, a 2 ] — > [— 1,1] 2 . The corresponding trajectory is a Lipschitz continuous map x(-) : 
[01,02] — > M such that x(t) = F(x(t)) + Ui(t)Gi(x(t)) + u 2 (t)G 2 (x(t)) for almost every 
t G [01,02]. Since the system is autonomous we can always assume that [01,02] = [0,T]. 

A solution to problem (P) is an optimal synthesis that is a collection {(xx{-),Ux(-)) defined 
on [0, Tx],x G M} of trajectory-control pairs such that x x (0) G M in , x x (T s ) = x, and 
x x (-) is time optimal. 

We use the following definition to describe different types of controls. 
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Definition 2. Let u(-) = (/Ui(-), v>2(-)) '■ [01,02] C [0,T] — > [— 1, l] 2 &e a control for the 
control system 01). 

• «(•) zs a frang control if for almost every t G [01,02], 



Similarly, u(-) is a unhang control (i = 1,2) if for almost every t G [01,02], Ui(t) = 



• A Ui(-)- switching (i = 1,2) is a time i G [01,02] such that for a sufficiently small e > 0, 
Ui(-) is a.e. equal to +1 on]i — e,i[ and a.e. equal to —1 on]i,i+e[ or vice-versa. A 
U1-U2- switching is a time i that is a u\- and a U2-switching. 

• If u A : [01,02] — > [— 1, l] 2 and Ub '■ [02,03] — > [— 1, l] 2 are controls, their concatenation 
Ub * Ua is the control 



The control u(-) is called bang-bang if it is a finite concatenation of bang arcs. Similarly 
one defines Ui-bang-bang controls. 



• A trajectory of (01) is a bang trajectory (resp. bang-bang trajectory) if it corresponds to 
a bang control (resp. bang-bang control). Similarly, one defines Ui-bang and Ui-bang- 
bang trajectories. 

Given two vector fields X and Y, consider the following function 



Notice that the definition of A(X,Y) depends on the choice of the coordinate system, but 
not the set Q{X, Y) that is the set of points where X and Y are parallel. For problem 



u(t) = ue {(-1, -l), (-1, l), (l, -l), (l, l)}. 



u G {±1}. 




A(X,Y)(x) :=det {X(x)),Y(x)), x G M, 



and the set of its zeros 



Q(X,Y) :={xeM s.t. A(X,Y)(x) = 0}. 
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(P), the sets Q(G 1 ,G 2 ), Q(Gx, [F + G 2 , G 1 }), Q(G h [F - G 2 , Gx]), Q(G 2 ,[F + G h G 2 }), 
Q(G 2 , [F - Gx, G 2 ]), and Q(F, G u G 2 ) : = Q(F, Gx) n Q(F, G 2 ) n Q(G U G 2 ) are fundamental 
to the construction of the optimal synthesis by applying PMP, as will be seen in the next 
paragraphs. In fact, assuming these sets to be embedded one-dimensional submanifolds of 
M, we have the following: 



u\-u 2 switchings can only occur on the set Q(Gx, G 2 ). See Lemma II. 3 



The support of Ux-u 2 singular trajectories (called totally singular trajectories in the 



following) is always contained in the set Q(Gx,G 2 ). See Lemma II. 6 



The support of Wi-singular trajectories (that are trajectories for which the Mi-switching 
function identically vanishes, and for which u\ can assume values different from ±1, see 
next section for detail) is always contained in the set Q(Gx, [F + G 2 , Gx]) U Q(Gx, [F — 



G 2 ,Gx])- A similar statement holds for -^-singular trajectories. See Lemma II. 7 



Under certain conditions, one proves that ux can switch only once on a connected 
component of M\(Q(F, G 1 )UQ(G 1 , G 2 )UQ(Gx, [F, ±G 2 , Gl]) U Q(G 2 , [F, ±G U G2])). 
A similar statement holds for u 2 . See Proposition |II.8 



For problem (P), Pontryagin Maximum Principle says the following: 

Corollary ILL Consider the control system subject to (HO). For every (x,\,u) E 
T*M x [-1,1] 2 , define 

H(x,\,u) := (\F(x)) + ux(X,Gx(x)) + u 2 (X,G 2 (x)) + A . 

If the pair (x(-),u(-)) : [0, T] — > M x [—1,1] X [—1,1] is time optimal, then there exist 
a never vanishing Lipschitz continuous covector A(-) : t 6 [0, T] \— > \{t) G T*, t sM and a 
constant Ao < such that for a.e. t e [0,T]: 

i) : x(t)=^(x(t),X(t),u(t)), 

ii) : X{t) = -^(x(t),X(t),u(t)) = -(X(t),(VF + Ul (t)VGx+u 2 (t)VG 2 )(x(t))), 
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iii) : H(x(t), X(t), u(t)) = 7iM(x(t),X(t)), where Hm(x,X) := max{"H(x, A, u) : u E 

[-M] 2 }, 

iv) : H M (x{t),\(t)) = 0, 

v) : (A(0), T x ( )Mi n ) = (transversality condition). 

Definition 3. The real-valued map H is called PMP-Hamiltonian . A trajectory x(-) (resp. 
a couple (x(-),X(-))) satisfying conditions i), ii), iii) and iv) is called an extremal (resp. an 
extremal pair). If (x(-), A(-)) satisfies i), ii), iii) and iv) with Ao = (resp. Ao < 0), then it 
is called an abnormal extremal pair (resp. a normal extremal pair). 

B. Switching Functions 

In this section we are interested in determining when controls switch from +1 to —1 or 
vice- versa and when they may assume values in ] — 1, +1[. Moreover we would like to predict 
which kind of switchings can happen, using properties of the vector fields F, G\ and G 2 - A 
key role is played by switching functions. 

Definition 4. (Switching Functions) Let (x(-), A(-)) be an extremal pair. The correspond- 
ing switching functions are defined as 4>i(t) := (X(t), Gi(x(t))), i — 1,2. For later use, we 
also define (f>o(t) := (X(t), F(x(t))) . 

The switching functions <f>\ and <p2 determine when the corresponding controls switch from 
+ 1 to —1 or vice-versa. In fact, from the maximization condition iii) of Corollary II. 1 one 
immediately gets: 

Lemma II. 2. Let (x(-), A(-)) be an extremal pair defined on [0,T] and(pi(-) the corresponding 
switching functions. If <fii(t) ^ for some t e]0,T[, then there exists e > such that x(-) 
corresponds a.e. to a constant control Ui = sign{<pi) on ]t — e,t + e[. Moreover, if <pi(-) has 
a zero at t, and if <fii{t) exists and is strictly larger than zero (resp. strictly smaller than 
zero) then there exists e > such that x(-) corresponds a.e. to constant control Ui = — 1 on 
]t — e,t[ and a.e. to a constant control Ui = +1 on ]t,t + e[ (resp. a.e. to a constant control 
Ui = +1 on }t — e, t[ and a.e. to a constant control Ui = —1 on }t, t + e[). 
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Notice that on every interval where &(■) has no zero (resp. finitely many zeros), the 
corresponding control is Mj-bang (resp. Uj-bang-bang) . Another direct consequence of the 
PMP is: 

Lemma II. 3. Let x(-) be an extremal trajectory defined on [01,02] and i G]ai,a2[ be an 
U1-U2- switching. Then x(t) G Q{G\,G 2 ). 

We are then interested in differentiating 0j. By a simple computation one gets: 

Lemma II. 4. Let (x(-), A(-)), defined on [0, T] be an extremal pair and <fii(-) the correspond- 
ing switching functions. Then it holds a.e. 

Mt) = (X(t),([F,G 1 ]+u 2 (t)[G 2 ,G 1 ])(x(t))), (5) 
Ht) = (A(t), ([F,G 2 ] +ui(t)[G l ,G 2 ](x(t))}, (6) 

where [•, •] denotes the Lie bracket of two vector fields. 



From Lemma II. 2 it follows that m can assume values different from ±1 on some interval 



[ax, a 2 ] only if the corresponding switching function vanishes identically on this interval. 



Remark 3. Lemma I 1. 4 asserts that if in a neighborhood of a u\-switching we have that u 2 
is a.e. equal to +1 or a.e. equal to —1, then in that neighborhood 0i(-) is a.e. a C 1 function. 
A similar statement holds for 4> 2 (-). 



C. Abnormal Extremals 

The following lemma is again direct consequence of the PMP. It characterizes some prop- 
erties of abnormal extremals. 

Lemma II. 5. Let (x(-),A(-)) be an abnormal extremal defined on [01,02]. We have: 

1. Ift is a U\-u 2 - switching, then x(t) G Q(F,Gi,G 2 ). 

2. If t is a u\-switching and u 2 is a.e. equal to +1 or a.e. equal to —1 in ]t — s,i + e[ for 
some e > 0, then x(t) G Q(F ±G 2 ,Gi). 

3. If i is a u 2 -switching and u\ is a.e. equal to +1 or a.e. equal to —1 in ]t — e,i + e[ for 
some e > 0, then x(t) G Q(F ± G%, G 2 ). 
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D. Singular trajectories 



Definition5. Consider an extremal trajectory x(-) defined on [01,02] • It is called u^- singular 
if the corresponding switching function </>«(■) vanishes identically on [01,02]- It is called 
totally singular if <j>i(-) and <j> 2 (-) both vanish identically on [01,02]. 

The following two lemmas are obtained immediately from the PMP. 

Lemma II. 6. Letx(-) be a totally singular trajectory on [01,02] C [0, T], then Supp(x(-)\[ ai) a 2 ]) C 
Q(G h G 2 ). 

Lemma II. 7. Let x(-) be a u\-singular trajectory on [01,02] C [0, T] and assume that u 2 is 
a.e. equal to +1 (resp. a.e. equal to —1) on [oi,a 2 ]. Then Supp(x(-)\[ aita2 ]) C Q(G\,[F ± 
G 2 , Gi\). Similar result holds true for u 2 -singular trajectory. 

E. Predicting switchings 

Definition 6. A point x G M is called a U\- super- ordinary point if 

x $ Q(F, G x ) U Q(G?i, G 2 ) U Q(G 1: [F ± G 2 , Gl]). 

On the set of u\- super- ordinary points we can define the functions ati(x), Pi{x), oji(x), ^i(x) 
as: 

[F, G x \{x) = a 1 {x)F{x) + f3 1 (x)G 1 (x), (7) 
[G 2 , Gi](x) = cu 1 (x)G 1 (x) + Ci(x)G 2 (x). (8) 

The following Lemma is not used in the rest of the paper, but is a step towards under- 
standing systems of the form (|4]). 

Proposition II. 8. Let Q C M be an open connected set composed of U\-super-ordinary 
points. Assume that for every x G Q we have ai(x) > and £i(x) — «i(x) > 0. (resp. 
ol\{x) < and C,i(x) — a\(x) < 0). Then all extremal trajectories x(-) : [01,02] — > fi, are u\- 
bang-bang with at most one —1 — > +1 u\-switching switching (resp. +1 — > —1 u\-switching) . 
A similar result holds for u 2 - switchings. 
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Proof. Let x(-) : ]ai,02[— > be an extremal trajectory and 0i(-) be the corresponding 
^-switching function. If 0i(-) has no zero, then x(-) is a Ui-bang and the conclusion follows. 
Let i be a zero of 0i(-)- The time i cannot be a zero of <f) 2 {-), otherwise we would have 
x(t) G Q(Gi,(j2)- From Remark [3] it follows that 0i(-) is a.e. C 1 in a neighborhood of t. 
Without loss of generality we assume that <pi(-) is C 1 in a neighborhood of t. Moreover i 
cannot be a zero of 0i(-) otherwise could not be a Wi-super-ordinary point (we would 
have x(t) G Q(Gi, [F ± G 2 , Gi])). Since in a neighborhood of i, u 2 is a.e. constantly equal 
to +1 or —1, we can assume u 2 constant in this neighborhood, and we have 

Mt) = (X(t),([F,G 1 ]+u 2 (t)[G 2 ,G 1 })(x(t))) 

= ax(X(t),F) +p 1 (\(t),G 1 ) +u 2 (t)u 1 (X(t),G 1 ) +u 2 (t)Z 1 (\(t),G 2 ) 
= axMt) + (Pi + Mtpi)Mt) + u 2 {t)£ 1( j) 2 {t). 

At time t, <f>\ = and if ai > and ^1 — ai > on Q, we have 

Mi) = ai(0 o (t) + u 2 (t)Mt)) + (6 - ocx)u 2 {t)Mt) > 0, 

where we used the following facts: i) from the maximization condition the quantity 
u 2 (i)4> 2 (i) > in a neighborhood of £; ii) Ao < implies that <fio(i) + u 2 (t)4> 2 (t) > —5 
for some arbitrary 5 > in a sufficiently small neighborhood of i (depending on 5). The 
case oil < and ^1 — ai < on Q is treated similarly. □ 



III. PROPERTIES OF EXTREMALS FOR SYSTEM (??) 

Based on the general results presented in Section [TTJ we derive properties of extremals 
for system ^ under assumption (Al), i.e., < a < n/A and < (3 < tt/4. We show 
in particular that starting from the north pole, only normal bang-bang trajectories can be 
optimal. All the results presented in this section are essentially based on the following lemma 
which characterizes the time evolution of switching functions corresponding to system (j3|. 

Lemma III.l. Let 4>$, <f>i, and <p 2 be the switching functions for system ((31). We have: 
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ft) 



P( Ul (t),u 2 (t)) 



, with 



P( Ul {t),U 2 (t)) 



( 

sin 2 a sin (3 cos /3 



cos a 



sin a sin f3 cos /3 



V 



cos a 



u 2 {t) 

Ui(t) 



(ii) On a bang-bang trajectory, we have 



cos a 
tan j3 



cos a 
tan j3 



u 2 (t) — cos a tan u\ (t) 
cos a tan /3 




\ 



0o(t) + |0i(t)| + |0 2 (t)| + A o = O. 



tan 2 a \sin 2 /3 cos 2 /3 



ifj.aia, w l thK:=^(^l+ m 



tan 2 a lsin 2 /3 cos 2 (3 



Proof, (i) is a consequence of Lemma II. 4 (ii) is a consequence of (iv) of Corollary II. 1 and 
Lemma II. 2 (iii) is based on (i) and the fact that 0o(O) = 0. □ 



A. Normal and abnormal bang-bang extremals 

Proposition III. 2. Under (Al), normal extremals for ^ have the following properties: 

(i) Let s > and s + t (t > 0) be two consecutive switching times. If <p 2 (s) = (resp. 
x (s) = 0), then ^ 0, + t) = and 2 (s + i) ^ ("resp. 2 (s) ^ ; 
2 ( s + t) = and +t)^0). 

(ii) The duration of the first bang-arc s satisfies [0,s max ] with 

sin 2 a cos 2 (3 



arccos 



arccos 



1 — sin 2 a cos 2 (3 
sin 2 a sin 2 (3 

1 — sin 2 a sin 2 /3 



corresponds to control (1,1) or (-1,-1), 
if it corresponds to control (1,-1) or (-1,1). 



(9) 
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(iii) The duration between two consecutive switchings is the same for all interior bang arcs 
(i.e., excluding the first and the last bang arcs). This duration depends only on the 
duration of the first bang arc. 



The proof of Proposition III. 2 is postponed to Appendix [Aj In the following, the duration 
between two consecutive switchings of interior bang arcs is denoted by v(s) with s the 



duration of the corresponding first bang arc. Point (iii) of Proposition III. 2 is illustrated in 
Fig. |3l and the explicit expression of v(s) is given in Appendix |Bj 




FIG. 3. Duration between interior switchings as a function of the first switching time s, for (3 = 7r/8 
with a = 0.25 (left) and a ~ ir/4 (right). The blue curve corresponds to initial control u = ±(1, 1) 
and the red to u = ±(1, —1). 



Corollary III. 3. If < a < it /A and (3 = 7r/4 ; we have 

(i) "5 max 



(a) v(s) 



arccos 



arccos 



sin 2 a 



1 + cos 2 a / 
d-A(s)-B(s)-C(s) 



, where A(s) := 8 cos a sin 2 a sin(s) ; B(s) 



e-A(s) + B(s) 

2 sin 2 2a cos(s) ; C(s) := 4 sin 4 a cos(2s) ; d := sin 2 2a, and e := 5 + 2 cos 2a + cos 4a. 



(iii) v(0) = v(s 



max/ «'iuux- 



(iv) All the switching points of the extremals having their first switching at s max are located 
on the great circles passing through N and containing x\— or X2—axis. 



Proof. Point (i) is a direct consequence of Point (ii) of Proposition III. 2 Points (ii)-(iv) are 
proved in Appendices [B] and [Cj □ 
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The next proposition concerns abnormal extremals. It holds without Assumption (Al). 



Proposition III. 4. There are no abnormal bang-bang trajectories starting from the north 
pole. 



Proof. Assume by contradiction that there exists an abnormal bang-bang trajectory starting 



from the north pole. Then, (ii) of Lemma III. 1 implies that 



Mt) + \Mt)\ + \Mt)\ = o- 



Since 0o(O) = 0, we have 0i(O) = 02(0) = 0. This contradicts the non triviality of the 
co- vector A. □ 



B. Singular trajectories 

The results presented in this section characterize singular trajectories of (|3]). They are 
consequences of Lemmas II. 6 , IL7[ and III.1[ The normalization used here for the co- vector 



A is given by A(0) = (cos#, sin#, 0) with 9 G [0, 2tt[. The corresponding initial conditions for 
the switching functions are: 

0o(O) = 0, (10) 
0i (0) = — sin a sin /3 sin 6* , (11) 
02(0) = — sin a cos (3 cos 9. (12) 



With this normalization, the constant K in (iii) of Lemma III. 1 is equal to cos 2 a. Moreover, 



it follows from (ii) of the same lemma that min^ Ao = — sin a. 

Proposition III. 5. The sets where the support of singular trajectories should belong to are 
characterized as follows. 

(i) The support of a totally singular trajectory must be contained in the equator Co of 
S 2 . The corresponding totally singular control satisfies U\ = and U2 = almost 
everywhere. 
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(ii) The support of a u\-singular (resp. U2-singular) trajectory must be contained in the 
set C\± := S 2 fl {± tan a cos (3 x 2 = x 3 } ( resp. C 2 ± '■— S 2 fl {± tan a sin (3 X\ = —x 3 } ). 
The corresponding U\-singular (resp. u 2 -singular) control satisfies a.e. U\ = and 
U2 = u, where u G {±1} (resp. u\=u and u 2 = a.e.j. 



Proof. For (i), applying Lemma II. 6 Gix(t) must be parallel to G 2 x(t). Therefore, x 3 (t) = 
on [a, b], i.e., a totally singular trajectory can only stay on the equator of S 2 . For (ii), assume 
for instance <f>\ — and 02 7^ on some interval [a, b]. Applying Lemmas II. 7 and III. 1 
G\x(t) is parallel to (G 2 — u 2 tan 2 a cos 2 (3 F)x(t), i.e., 



/ o \ 



( u 2 tan a cos f3 x 2 — x^ 



A 



V 



-u 2 tan a cos (3 x± 



(m 2 tan a cos a x 2 — x 3 )x = 0. 



-x 3 
\ X2 / 

Therefore, a Mi-singular trajectory must stay in the set C\± := S 2 fl {± tan a cos (3 x 2 = x 3 }. 
The proof for ^-singular trajectory is similar. □ 




FIG. 4. Co, C\±, and C 2 ± are represented respectively by green, red, and blue lines for a = ir/5 
and f3 = tt/6. 
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Corollary III. 6. Assuming (Al), a normal bang-bang extremal starting from the north pole 
cannot connect to a totally singular arc. 

Proof. Assume by contradiction that a normal bang-bang extremal is connected to a totally 



cos 2 a, 



singular extremal at t c . Then by Lemma III. 1 , we have 

0o (*c) = -Ao, 4>l{t c ) = 
implying that Aq = cos 2 a. However, we already now that 

max A 2 , = sin 2 a < cos 2 a, if a < ir/4, 



□ 



which yields a contradiction. 

Corollary III. 7. If a < tt/A, a normal bang-bang extremal starting from the north pole 
cannot connect to a partially singular arc. 

Proof. Assume for instance <p2(t c ) = 4>2(t c ) = and <j)i(t c ) ^ 0. Note that if <pi(t c ) = 0, then 



we proceed as in the proof of Corollary III. 6 to achieve a contradiction. From (ii) and (hi) 
of Lemma [HI. 1 , we obtain 



0o(t c ) + |0i(* c )| = -A o , 



<f>p(tc 

cos 2 a ' sin 2 a sin 2 6 



+ 



4>m 



which imply 



00 (t c ) 



\MQ\ 



— Ao (sin 2 a sin 2 (5 — cos 2 a) =F 

sin 2 a sin 2 (3 + cos 2 a 
— A sin 2 a sin 2 (5 ± \/A 
sin 2 a sin 2 (5 + cos 2 a 



where A := 1 — An — sin 2 a cos 2 (3. From (i) of Lemma III. 1 , we also have 



Mtc) 



sin 2 a sin (3 cos (3 



Ui(f>o(t c ) - 



cos a 



tan/5 



01 (tc) 



cos a 

sin 2 a cos 2 B + cos 2 a „ r— 
=F — n sin/3MiVA = 0, 



cos a cos /3 

which implies A = 0. However, if a < 7r/4, then Ag < cos 2 a and 

A > sin 2 a sin 2 /3 > 0, 



which yields a contradiction. 



□ 
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IV. OPTIMAL SYNTHESIS FOR SMALL a AND 8 = - 

4 

Assuming (A2), i.e., a small and (3 = it /4, we prove in this section that the extremal 
trajectories issued from the north pole are globally optimal until they reach a neighborhood 
of the south pole. In the following we use S 2 \ 0(a) to denote the sphere S 2 minus a 
neighborhood of order a of the south pole. 

In general, proving global optimality of solutions of the PMP is not an easy task since 
one has to compare for each final point all extremals reaching that point. In our case, we 
get the result by a set of arguments similar to those used for the problem with one bounded 
control, see 0HI)]- We only give a sketch of these arguments (to avoid lengthy computations 
similar to those made in [7J [15]) except for one crucial proposition that is proved in detail. 
These arguments are described in the following steps. 

STEP 1 We consider all extremals starting from the north pole. They are divided into 4 families 
depending on the value taken by the controls at the beginning, namely (1, 1), (1, —1), 
(-1, -1), (-1, 1). Let A" sign ( Ul ) sign ( U2 ) := F + u i G 1 + u 2 G 2 . Then, the first family of 
extremals has the form: 

E(t, s, a) := m(t, s, a)M n {s, a)e sX++ N, s G [0, s max ], (13) 

where n is an integer, M(s, a) is defined by 

M(s, a) := e <»)x ++e v(,)x- +e v{.)x— e v(s)x+- i ^ 

and m(t, s, a) has one of the following forms, 



m(t, s,a) = < 



e (t-n(n, a ))A-+- j n ( n) s ) ._ 4 nw ( s ) _|_ S) 

e (t-r 2 (n,s))X— e v(s)X+.^ g J ;= ^ + + ^ 

e (t-T 3 (n,s))x„ +e v(s)X-_ e v(s)x+_ ^ r 3 (n,s) := (4n + 2)v(s) + s, 

e (t-T 4 (n,s))X ++e v(s)X- + e v(s)X— e v(s)X + _^ T ^ ^ ._ ^ + 3^^) + ^ 



with < t — Ti(n,s) < v(s). Note that the integer n and the function m(t,s,a) are 
determined by the target. All the three other families can be defined in a similar 
manner. The extremal trajectories having their first switchings at s max are called 
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"boundary-trajectories" of the family. Each extremal of the family switches a certain 
number of times before reaching the south pole, and all the switching points of the 
different extremals form smooth curves called switching curves. See Figure |6j 

Definition 7. A (geometric) smooth curve C is called a switching curve if each point 
of C is a switching point. 

By construction, the following curves Ck{s, a) defined by induction are switching curves 
of the first family: 



d(s,a) = e sA ++iV, C k (s,a) = M(s,a)C k ^(s,a), 



(15) 



with k > 1, s G [s,s max ], and M(s, a) defined by (14). 




switching curve 
refracting the trajectories 



switching curve 
reflecting the rajectories 



switching curve 
reflecting and refracting 
the trajectories 



FIG. 5. Local optimality of switching curves 



STEP 2 Each switching curve can "refract" or "reflect" the extremals (see Fig. [5]). The main 
argument of the proof is that up to a neighborhood of the south pole, all the switching 
curves are "locally" optimal, i.e., they always "refract" extremals. 

Definition 8. Let C be a switching curve. Assume that extremal trajectories switch 
on C from a smooth vector field Y\ to another smooth vector field Y 2 . Let C(s) be a 
smooth parameterization of C with s G Domain(C) . We say that C is locally optimal 
if, for every s G Domain{C) and for every pair (ci,c 2 ) such that CiC 2 > 0, we have 
d s C(s)^c 1 Y 1 (C(s)) + c 2 Y 2 (C(s)). 
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Proposition IV. 1. Let a be small enough and consider the set of extremals issued 
from the north pole before they reach a neighborhood of the south pole. Then all the 
switching curves formed by this extremal flow are locally optimal. 

The proof of this proposition is given in detail at the end of the section. It is clear 
from numerical simulations that the neighborhood of the south pole where extremal 
flow loses optimality is approximately a disk of radius 3a. 




*■ = trajectory 
\ = switching curve 



FIG. 6. The four-snake structure. 



STEP 3 Proposition IV. 1 has two main byproducts (which are not completely obvious, but can 



be proved as in H3] ) . 

— The four families of extremals defined in STEP 1 are well organized in a structure 
of four snakes (see Fig. [6] and Fig. [T]), in the sense that the four families do not 
intersect until they reach a neighborhood of the south pole. 
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— In each snake, trajectories do not intersect each other until they reach a neigh- 
borhood of the south pole. 

As a consequence, each point of S 2 \ 0(a) is reached by one and only one extremal 
issued from the north pole before reaching a neighborhood of the south pole. By 
construction, these trajectories are optimal. 

The four snakes intersect in a neighborhood of order a of the south pole. Hence the proof 
fails in that region. One can also see that in a neighborhood of the south pole there exist 
non locally optimal switching curves. An analysis on how the trajectories lose optimality in 
a neighborhood of the south pole is very complicated and out of the purpose of this paper, 
but it can be pursued as in [15]. The extremal front, defined as the set of the endpoints 
of extremal trajectories at time t, is homeomorphic to a circle up to a neighborhood of the 
south pole. In a neighborhood of the south pole it develops singularities (cusps and self- 
intersections) showing the presence of a cut locus (locus at which trajectories lose optimality). 
See Figure [7j 




FIG. 7. Top view of the extremal front, far from the south pole (left) and close to the south pole 
(right). 



Proof of Proposition IV. 1 By (15), we obtain 

C k+1 (s,a) := M(s,a)C k (s,a) = M k (s, a)d(s, a). 
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(16) 



The curves Ck(s,ct) correspond to a switching from (+, +) to (+,—). Assume that there 
exist two real numbers C\ and c 2 such that Cic 2 > 0, and 

d s Ck+i = (c\X ++ + C2-X" + _)Cfc+i. 



Therefore, 



(ciX ++ + c 2 X + _)M fc Ci = M k d s d + M^dsMM^d. 



i=i 



By Taylor expansion, we also have 



(17) 



Ci(s,a) 



9 s Ci(s,a) 



M(s,cc) 



— sins — cos sj ct 

Ai ^ • ^ 

--— (1 + sins — cos s) a 



\ 1 

— (— cos s + sin s) a 

A ± ■ \ 

- —— (cos s + sin s) a 



+ 0(a 2 ), 



I 







+ 0(a 2 ) 



/ 



\ 



+ 0(a 2 ) = i?((9) + 0(a 2 



<9 s M(s,a) = 



1 
1 -A^2a 
[0 4y/2a 1 J 

4(sins — coss) a 2 0^ 

—4 (sins — cos s) a 2 





+ 0(a 3 ) = 0(a 2 ), 



with R(6) :-- 



1 
cos 6 — sin 9 
sin# cos# 



\ 



and 9 := 4^2a. 
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By (17), we obtain 



(Cl + c 2 



sin kO^ 


V I 



( v% 

d s fi(s) + (ci + c 2 ) cos k6 f 2 (s) + (ci - c 2 )— cos k8 

cos k9 d s f 2 (s) - (ci +c 2 )[/i(s) — cos k6] 

sin k6 d s f 2 + (ci + c 2 )— - sin 



V 



a + 0(a 2 ), 



/ 



y/2 a/2 
with /i(s) := ~2~(1 — sins — coss) and / 2 (s) := — — (1 + sins — coss). Eq. (18) implies 

that 

sin k6 = sin 4^2A;a = 0(a). (19) 

Clearly this condition can be satisfied only in a neighborhood of order a of the north pole 
or of the south pole. Direct computation shows that it is not satisfied in a neighborhood of 
the north pole. It follows that the switching curves Ck are locally optimal until intersecting 
a neighborhood of order a of the south pole. All the other cases can be treated in a similar 
manner. □ 

V. SIMPLE SUBOPTIMAL CONTROLS AND COMPARISON WITH OTHER 
STRATEGIES 

A. Two simple suboptimal strategies realizing complete spin flip 



Based on the optimal synthesis described in Sec. IV and the computational lemmas 
gathered in Appendix [Cj we present in this section two simple suboptimal strategies for the 
case where the two controls have the same bound M and the ratio M/E is small, i.e., a is 
small and (5 = tt/A. This case is the most relevant one for NMR applications. 



1. A really simple suboptimal strategy 



We first present the strategy for the normalized system ^ with k = 1. For small a, we 



obtain from Lemma C.4 that v(s) ~ tt/2. Consider the following sequence of controls 
(S 1 ) : (+1,-1) -> (-1,-1) (-1,+1) -> (+1,+1), 
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where each combination of (111,112) lasts for a duration of tt/2. By Lemma C.5 (see also the 



proof of Proposition IV.l), the action the sequence (Si) produces approximately a rotation 
around x\— axis of angle 4y / 2<x Let n := [ — -j=- ] , where \K~\ denotes the smallest integer 



7T 



4\/2a 

< a < 



71 



-. Then, starting from the north 



not less than K . It is clear that ,_ 

Ay/2n Ay/2(n - 1) 

pole, applying the sequence Si for n times steers the system close to the south pole (see Fig. 



8|), and the error is of the order of a. 




2 4 6 8 10 12 

Time 



FIG. 8. On the left u\ (top) and 112 (bottom) over two periods. On the right, suboptimal trajectories 
for a = 0.01 and k = 1. 



We now take into account the time normalization constant k which is approximately equal 
to 2E in the case of small a. The suboptimal controls corresponding to the sequence Si are 
periodic rectangular signals of period tt/E, and the total transfer time is equal to nir/E. The 
advantage of this strategy is that it only requires the knowledge of the Larmor frequency 
of the system and the bound on the control fields; it does not necessitate any computation. 
Note also that the three other sequences obtained from cyclic permutations of (Si) are also 
suboptimal in a similar manner. Finally, this strategy is suboptimal in the sense that the 
south pole is not exactly reached. 
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2. A more accurate suboptimal strategy 



We investigate more carefully properties of the sequence (Si) to derive another suboptimal 
strategy that steers (|3| from the north pole to the south pole exactly. The new strategy is 



based on the following proposition, which is a direct consequence of Corollaries |C.2| and |C.3 

-2\/2 sin a cos a 



Proposition V.l. Let 



la 



arcsm 



1 + cos 2 a 



, and M(s, a) be defined by Eq. 



(14). For any interger n and any a < |, we have 
( 



M?(Q,a)N 







\ 



sin(4n#(a;)) 
^cos(4n6'(a))y 



/ 



M?(s t 



a e 



-N 



\ 



sin((4n + l)O(a)) 
\^cos((4n + l)9(a)) J 



(20) 



In other words, Proposition V.l states that for any a less than ir/A the switching curve's 
two endpoints stay on the great circle orthogonal to the X\— axis under the action of the 
sequence (Si). We now construct a suboptimal strategy from the formula for M"(0, oc)N. 

71 

Let n := ] . It is easy to check that 9(a) is a monotonically increasing function of a 

for a G [0,7r/4]. Then, there exists a unique a < a such that An6(a) = it. In other words, 
we have 



M?(0,a)N 



/ N 



v-v 



(21) 



Let 7 := 



sequence 



sin a 
sin a 



< 1. Then, the action of Mi(0,a) is realized by applying the following 



(S 2 



(+7, -7) -> (-7, -7) -> (-7, +7) -> (+7, +7), 



where each combination of controls lasts for a duration equal to arccos 



sin 2 a 



1 + cos 2 a 



which 



is approximately tt/2. This strategy is suboptimal in the sense that the transfer time is not 
optimal, because the controls do not satisfy the Pontryagin Maximum Principle (the controls 
are bang-bang, but do not saturate the bounds ±1). The advantage of this strategy is that 
the south pole is exactly achieved for any a less than tt/4, and the transfer time is close 
to the optimal one if a is small. See Fig. M for a comparison with the strategy (Si) for 
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different values of a. Three other similar suboptimal sequences can also be built from cyclic 
permutations of (S 2 ). 




FIG. 9. Trajectories of X3 corresponding to strategy S\ (blue line) and strategy 52 (red line) for 
a = 0.1 (left) and a = 0.01 (right). 



B. Comparison of suboptimal controls with optimal control bounded on the circle 

In this section we compare the times needed to steer system ^ from the north pole to 
the south pole for two different bounds on the controls: 

{S): \Qi{t)\ < M, i = 1,2, (CO : ^ Ql(t) + Q 2 2 {t) < M. 

We are only interested in the case where M is much smaller than E, i.e., a is small. It was 
shown in [5], [6] that the following control is time optimal for the control system with control 
bound of type (C): 

nUt) = Msm(tu r t + 0) , 

(22) 

Q 2 (t) = Mcos (u r t + </)), 

where u r = 2E, and that the optimal time to steer from the north pole to the south pole is 
given by 

Tc(M) = ^. (23) 
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We now estimate T$(M), the optimal transfer time for (S). Recall that if Ts(a) is the optimal 
time for the normalized system ^ with k — 1, then we have Tg(M) = Te(ck) / (2a/ E 2 + 2M 2 ). 



Using for instance the first suboptimal strategy presented in Sec. V A (the second suboptimal 
strategy gives the same result up to an error of order a), we obtain 

■2 



Tv(a) = 2tt 



7T 



+ o(a) 



7T 



+ o(a). 



AV2aJ 2V2a 
We also have 2a/ E 2 + 2M 2 = 2E + o(a). Therefore, T$(M) is approximately given by 



(24) 



T S {M) 



7T 



7T 



4A/2aE 8M 



Eqs. (23) and (25) imply that 



T S {M) 7T 



0.78. 



(25) 



(26) 



T C (M) 4 

In other words, there is an improvement of 22% when using suboptimal controls for problem 
(S) compared to the optimal ones for problem (C). 




FIG. 10. Time evolution of (u\,U2) in the limit a — > for the control system on the square (red) 
and the circle (black). 



Appendix A: Proof of Proposition III. 2 



Recall that the normalization for the co-vector A is given by A(0) = (cos 9, sin 6, 0), and 



the corresponding initial conditions for the switching functions are given by Eqs. (11) and 
(12). It follows from (ii) of Lemma III. 1 that mine Aq : 



sin a. 
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Proof of (i). By (ii) and (iii) of Lemma III. 1 we have 



0o(t) + |0i(t)| + |0 2 (t)| + A o = O 



cos 2 a sin 2 a sin 2 j3 sin 2 a cos 2 j3 



+ 



+ 



2 (t) 



(Al) 
(A2) 



We call the sets of (0o, 0i,02) that satisfy equations (Al) and (A2), Su and S a d respec- 



tively. By (i) of Lemma III. 1 the solution of the adjoint system is defined and unique in 
[s, s + 1]. This solution must always lie in the intersection S-n R Sad, which is visualized in 



Fig. [11} The surface defined by S% is a union of four quarter-planes and may as such be 
spanned by rays starting at their common intersection, 



-An 



0. 



0. 



(A3) 



FIG. 11. Intersection of the surfaces S-u and S a( i- 



Assume by contradiction that <p2 (s) = 02 (s + t) = 0. This implies that there exists a 



ray in starting from (A3) that intersects S a d more than once, since the solution must 



return to the (0o ; </>i)-plane without crossing itself. However, if a < 7r/4, (A3) belongs to 



the interior of S a d, which is a strictly convex set. Thus any ray starting from (A3) can only 



intersect S a d once. Thus 02 can never switch two times consecutively. The same reasoning 
holds for 0i. □ 
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Proof of (ii). Assume that the initial control is (1,1), i.e. 9 G [tt,3tt/2]. By (i) of Lemma 



III. 1 with the initial condition given by A(0), for a given 9 G [k, 37r/2], s is the first switching 



time if and only if f(8, s) = with 



. . (2 (1 -sin 2 a cos 2 /?) . 
f(9, s) := 2 cos @ — s i n 2/3 sin 6 1 cos s 

1 sin a ' 

2 COS OL 

sin # sin s + (1 + cos 2(3) cos # + sin 2/3 sin 



sin 2a 

sin 2 a cos 2 /3 

1 — sin 2 a cos 2 (3 
problem 



It follows directly that coss(7r) = — . Consider the following initial value 



s 



s(tt) = arccos 



d s f(e, s 

sin 2 a cos 2 (3 



a«_a« x ( A4 ) 



1 — sin a cos 2 /3 



where s' denotes the derivative of s with respect to 9. It is easy to show that (A4) has a 
unique solution for 9 G [it, n + e[ with e > small enough. 



The first step consists of showing global existence of the solution to Eq. (A4) over the 

interval [tt,3/2tt]. We have 

(2 (l - sin 2 a cos 2 /?) \ 2 cos a 

defp(6, s) — — 5 sin 9 + sin 2/3 cos 9 cos s , — cos # sin s 

\ sin a / sin a 

- (1 + cos 2/3) sin 9 + sin 2/3 cos 9, 

(2 (l - sin 2 a cos 2 /3) \ 2 cos a 

o s fp(6, s) — — g cos 6* — sin 2/3 sin # sin s ^ — sin 9 cos s. 

\ sin a / sin a 

It is clear that there exists a constant K > such that |<96i/(#, s)| < ii". We need a uniform 
estimate for l/d s f(9,s). Let 

2(1 — sin 2 a cos 2 B) , N 2 cos a 

a(0) := , cos - sin 2/3 sin 0, b{9) := 5— sin 9, 

sin a; sin a 

c(0) := (1 + cos 2/3) cos 9 + sin 2/3 sin 9. 

We have that /(#, s) = a(9) cos s — b(9) sin s + c(9). It follows that 

/(M • • , c(G) 

cos 7 cos s — sin 7 sin s 



y/a 2 (6) + lP(e) y/a*(6) + b 2 (6) 

c(9) 

= costs + 7H . , 

y/a*(6) + b*(6) 
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where cos 7 :- 



a{9) 
^a 2 {9) + b 2 {9) 



and sin 7 :- 



^a 2 (9) + b 2 (9)' 



Therefore, \d s f p (9,s 



a 2 (9)+b 2 (9) |sin(s + 7 )|. 



Since f(9,s) = 0, we also have that 



sin 



> + 7) I = a/1 - cos 2 (s + 7) = a 1 - 



c 2 {9) 



a 2 (9) + b 2 (9) 



which implies 



\d s f p (9,s)\ = y/a 2 (9) + b 2 (9)-c 2 (9). 

We simplify the expression within the square root, and obtain 

2 



a 2 (9) + b 2 (9)-c 2 {9) 
We deduce that 



sin 2 a 



(2 cot 2 a - 2 cos 2/3 cos 2 9 - sin 2(3 sin 29) . 



\d s f] 



S)\ > 



2-y/cot 2 a — cos 2 (3 



and it follows that 



deft 



d s f P { 



< 



sin a 



K sin a 



2\/cot 2 a — cos 2 (3 
which implies that the solution of (|A4[) is globally defined on [tt, 3tt/2]. 



(A5) 



(A6) 



Let s(-) be the solution to Eq. (A4). The second step consists of showing that s is 



decreasing on [tt, 3ir/2]. We deduce from (A5) that d s f(9, s) does not change sign on [ir, 3ir/2] 
and it is easy to show that it is positive. We now show that dgf(9, s) does not change sign 
neither. By contradiction, assume that dgf(9, s) = for some s. Together with the fact that 
f(9, s) = 0, we obtain 



with 



-c\ sin 9 — C4 cos 9 —C2Cos9 





' cos s 1 


H 




I sin s 1 





ci :-- 
c 2 := 



2(1 



sin a cos* 



0) 



2 cos a 
sin 2 a 



sin 2 a 



c 3 := 
c 4 : = 



— C3 cos 9 — C4sin# 
C3 sin 9 — C4 cos ^ 

1 + cos 2/3, 
sin 2/3, 



(A7) 
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which implies 



cos s 



sin s 



sin 2 a cos 2 (3 
1 — sin 2 a cos 2 (3 ' 
sin 2 a sin (3 cos /? 



(1 



sin a cos^ p cos a 



Therefore, 



cos 2 s + sin 2 s 



sin 2 a cos 2 (3 



(1- 



sin a cos^ p cos z a 



(A8 



However, it is easy to show that ( A8 ) is less than one if a < 7r/4, thus yielding a contradiction. 



Therefore def(8, s) does not vanish and s' has a constant sign on [it, 37r/2]. Using 

— sin 2 a cos 2 (3 



COS S 7T 



sins 7T 



1 — 2 sin a cos 2 (3 



1 — sin 2 a cos 2 j3 ' 1 — sin 2 a cos 2 /3 ' 

it is easy to show that dgf(ir, s(tt)) > 0. Therefore s'(tt) < 0, which implies that s' < on 
[7r, 37r/2l. We conclude that max s(9) = sOk). The other cases with initial control equal 

6»e[7r,37r/2] 

to (—1,-1), (—1,1), or (1,-1) are similar. □ 



Proof of (hi). Assume for instance </> 2 (s) — (pi( s + ti) = ^(s + 1\ + t 2 ) = 0. Let X Ul , 
F + MiGi + w 2 G 2 . Then, we have 

\(s + t x ) exp(tiX UlU2 )G 2 exp(-tiX UlU2 )2;(s + t x ) 
= A(s + *i)G r i«(s + t 1 ) 

= A(s + h) exp(-t 2 X_„ lU2 )G 2 exp(t 2 A_ UlU2 )2;(s + t x ) 
= 0. 



(A9) 



Recall that the Lie algebra (so(3),[, ]) is isomorphic to the Lie algebra (IR 3 ,A), where A 
denotes the vector product in M 3 , and we use the following isomorphism: 



/ 



i : 



-c b 
c —a 
-b a 



6 

V7 



Then, (A9) is equivalent to 



det(i(Gi),e tlA "i"n(G 2 ),e 



*i(G 2 )) = 0, 
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which may be simplified to the following equation 



sm I \ D{tx,t 2 )=0, 



(A10) 



where 



D(t 1 ,t 2 ) ■= 2sin(|) (4 cos a cos (3 cos (|) +sin(|) ( (cos 2a + 3) sin (3 - 2 sin 2 a sin 3/3)) 
+ 8cos (y) (cos a cos (3 sin (|^) + sin f3 cos (tt)) . 

It is easy to check that D(ti,t 2 ) ^ if (£i,t 2 ) £ [0,tt] 2 , which implies that t x = t 2 . Lemma 



A.l guarantees that ti and £ 2 indeed belong to [0, 7r]. Therefore, the duration between any 
two switchings is the same. □ 

Lemma A.l. Let s, s + t\ and s + t\ + t 2 be three consecutive switching times. Then 
(h,t 2 ) G [0,tt] 2 . 



Proof of Lemma A.l . To simplify the notation, set a\ := cos a cot /3, a 2 := cos a tan /3, and 
a 3 := cos a tan 2 a sin (3 cos /?. We have 



P(ux,u 2 ) 



^ M 2 ai — Uia^ 
-u 2 a 3 a 2 
Wia 3 — ai 



and aia 2 + aia 3 + a 2 a3 = 1. 



Let 



u 2 a 3 
u\a 2 

u 2 ai 
V 1 / 



, v 2 (u 1 ,u 2 ) ■-- 



/ u 2 at \ 

ai + a 3 
u\u 2 a 3 

a± + a 3 

1 



V 



ai + a 3 
1 



It is straightforward that Pv\ = 0, Pv 2 = t> 3 , and Pv$ = —v 2 . 

Let Q(vi,v 2 ) := [f i(«i, u 2 ), v 2 {u\, u 2 ), v s (ui, u 2 )]. We have Q~ 1 PQ 



ai + a 3 




A) \ 
0-1 

v° 1 °/ 



with 



34 



/ 



Q- 1 



u 2 a\a 3 U1U2CL1CI3 a x a 3 
-U2CL1CL3 —U\U 2 aia 3 02(01 + 03 
—U1CL3 ai 

/ 

Q(ui,u 2 ) 



\ 



Therefore, 



1 

cos t — sin t 
sint cost 



Q (ui,u 2 ) 



(All) 



k{s + t) 

\M s + t)J 

where (ui, u 2 ) are the controls used in the interval [s, s + t}. 

Let s and s + t be two consecutive switching times. We show that t G [0, n]. Without loss 
of generality, we can assume that <f> 2 {s) = 0. Then, we have 0i(s + 1) =0. We also know 



that 4>o(s) = — Ao — U\(j)\{s). Therefore, (All) implies 

4>i (s + t) = (pcos a + Aq cos a sin 2 a sin 2 (3) cost + U\U 2 sin (3 cos (3 sin 2 a(p + Aq) sin t 



— A cos a sin 2 a sin 2 /3 = 0, 



where p := Ui0i(s) 



>i(s)| is the positive solution of 
1 



1 + 



tan a sin J3 



)p 2 + 2\ p + Aq - cos 2 a = 0. 



Note that (A12) has exactly one positive solution if a < tt/4. Therefore, we obtain 

-Aq sin 2 a sin 2 j3 + cos a sin a sin /3\/A 



P 



sin 2 a sin 2 j3 + cos 2 a; 



(A12) 



(A13) 



(A14) 



with A := 1 — A — sin a cos (3. Substituting (A14) in (A12), we obtain after simplification 



where 



ci :-- 



c\ cos t + c 2 sin t + C3 = 0, 



Aq sin (3 sin 2 a cos 2 (3, 



(A15) 



cosav A 



sma 



c 2 := M1M2 cos /3(sina sin/3v A + Ao cos a), 
c 3 := — A sin/3(sin 2 (3 + cos 2 a cos 2 /3). 



Eq. (A15) is equivalent to 



(c 3 - c x ) tan 2 - + 2c 2 tan - + (c 3 + c x ) = 0, 



(A16) 
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and we only need to show that (A16) has positive solution. Prove that - — - < 0, which 



c 3 - Ci 



implies that (A16) has exactly one positive solution. As Aq < 0, it is clear that C3 + c% > 0. 



Therefore, we only need to check the sign of C3 — c\. 

A sin a sin (3(1 — 2 sin 2 a cos 2 (3) + cos ay/A 



Let A 







c 3 - Ci 



A > 0. We have 



sma 



cos ay 1 — Aq — sin 2 a cos 2 (3 — Ao sin a sin (3(1 — 2 sin 2 a cos 2 /3) 
cos 2 a(l — Aq — sin 2 a cos 2 (3) — A 2 sin 2 a sin 2 — 2 sin 2 a cos 2 (3) 2 

cos oja/I — Aq — sin 2 a cos 2 /3 + A sin a sin — 2 sin 2 a cos 2 /3) 
cos 2 a(l — Aq — sin 2 a cos 2 /3) — Aq sin 2 a sin 2 — 2 sin 2 a cos 2 /3) 



cos a 



\J\ — Aq — sin 2 a cos 2 /3 + A sin a sin — 2 sin 2 a cos 2 f3) 



cos 2 ck(1 — sin 2 a cos 2 /?) — Aq(cos 2 a + sin 2 a sin 2 (3 — 2 sin 4 a sin 2 /3 cos 2 f3) 



> 



cos a a/1 — Aq — sin 2 a cos 2 /3 + Ao sin a sin — 2 sin 2 a cos 2 /3) 
Using the fact that, for a < ir/4, cos 2 a > Aq, and 

(1 — sin 2 a cos 2 /?) — (cos 2 a + sin 2 a sin 2 (3 — 2 sin 4 a sin 2 /3 cos 2 (3) 
= 2 sin 4 a sin 2 /3 cos 2 /3 > 0, 



(A17) 



we conclude that (A17) is positive, which implies that C3 — c\ < 0. 



□ 



Appendix B: Expression of v(s) 

We establish in this section the expressions for the duration between switchings on a 
normal extremal as a function of the the first switching time s. For convenience we set c a 
and s a to denote respectively cos a and sin a. 



Lemma B.l. // the extremal trajectory is starting from the north pole with control (1,1) 
or (—1,-1) (resp. with control (1,-1) or (—1,1)^ then v(s) = Vi(s) (resp. v(s) = V2{s)), 
where 

' 'A i {s)+B i {s) y /Q(s)TE~ i 



vas) = arccos 



DAs) + F l 



1,2, 
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with 

Ai(s) = 4^^2/3 + 8sl a sp(c a cp sin(s) + sp cos(s)) + 2s£ ((3 + c 2a )s|g cos(2s) + 2c a s 4 p sin(2s)) . 

A 2 (s) = 4:S 6 a sjp + Ss 2 2a cp (c a sp sin(s) + cp cos(s)) + 2s 4 a ((3 + c 2a )s|g cos(2s) - 2c Q s 4/3 sin(2s)) 

Si(s) = V2s 2 a cp (cp sin(s) + c a sp(cos(s) - 1)) , 

.62(5) = V2s 2 a sp (sp sin(s) + c a cp(cos(s) - 1)) , 

Ci(s) = 256s 2 a c a cp (3 + c 2a + 2s 2 a c 2 p) (c a cpCOs(s) - spsm(s)) 

+ 64^4 ((( 3 + c 2a)c 2 /3 - 4:S 2 a ) cos(2s) - 4c a s 2/ 3sin(2s)) , 
C 2 (s) = 256s 2 a c a sp (3 + c 2a - 2s 2 a c 2 p) (c a spCOs(s) - cpsin(s)) 

- 64^4 (((3 + c 2a ) c 2 p + As 2 a ) cos(2s) + Ac a s 2 p sin(2s)) , 
Di{s) = 16s 2 a c a cp (3 + c 2a + 2s 2 a c 2 p) (spsm(s) - c a cpCOs(s)) , 
D 2 (s) = 16s 2 a c a sp (3 + c 2a - 2s 2 a c 2 p) (c^sii^s) - c a sp cos(s)) , 

£1 = 234 + 3844c 2 ^ - 16s^c 4/3 (l + 3c 2a ) + 205c 2a + 70c 4a + 3c 6a , 

E 2 = 234 - 3844c 2/3 - 16s* 0^(1 + 3c 2a ) + 205c 2q + 70c 4a + 3c 6a , 

Fx = -17 - 16s^c 2 ^ + c 2q ( 4s*c 4/9 



F 2 = -17 + 16s*c 2/ 9 + c 2a 4s* C4/J 
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1 




— 5c 4q , — 




39\ 


— 5c 4a — 


1 




4 C6a ' 



Proof. Similar to the proof of (iii) of Proposition III. 2, assume 

(j) 2 {s) = <p 1 (s + v) =0, 
where s is the first switching time. Therefore, 

X(s)G 2 x(s) = A(s) exp(— vX + J)Gi exp(vX + _)x(s) = 0, 

which implies that 

g(s,v) := det(z(G 2 ), i(exp(— t>X + „)Gi exp(t>X + „)), exp(sX ++ )a;(0)) 

= 0. (Bl) 

After simplification, we obtain 

a(s) cosv + b(s) smv + c(s) = 0, (B2) 
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with 



a(s) := — cos 2 (3 cos s + cos a sin (5 cos (5 sin s — cot 2 a, 

b(s) := — cos a sin (3 cos /3 cos s — cos 2 /3 sin s + cos a sin (3 cos /3, 

c(s) := — sin 2 j3 cos s — cos a sin /3 cos sin s. 



The result follows. 



□ 



Lemma B.2. For the special case of (3 = n/4, we obtain a simpler expression: 

~d-A(s)-B(s)-C(s) 



v (s = arccos 



e-A(s) + B(s) 



where 



A(s) = 8c a s 2 a sin(s), d = s| a , 

5(s) = 2s 2 a cos(s), e = 5 + 2c 2a + c 4a , 

C{s) = 4s^cos(2s). 



(B3) 



From Lemma B.2 and (ii) of Proposition III. 2 we deduce the following corollary. 



7T 



Corollary B.3. For f3 = —, we have v(0) = v(s mSLX ) = s max . 



Proof. From (ii) of Proposition III. 2, we have 



sin 2 a 



cos s r 



1 + cos 2 a ' 



smsr 



2 cos a 
1 + cos 2 a 



(B4) 



Substituting (B4) into (B.2), it is easy to check that cos(i>(s max )) = cos(s max ). In a similar 
manner, we check that v(0) = s max . □ 



Appendix C: Some computational lemmas for the case j3 = tt/4 



The following result is a consequence of Corollary |B.3 It can be checked by direct 
computation. 
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, -2y/2 sin a cos a 
Lemma CI. Let 9(a) := arcsm | — ; ) . lhen, we have 



3 «(Smu)X++ 



0-10 
cos 0(a) sin0(a 
y— sin8(a) cos0(a)y 

'0-1 ^ 

cos 0(a) — sin 0(a) 
, sin 0(a) cos 9(a) j 

sin a cos a 



1 + cos 2 a 
\ 



max )x+. 



,»(s milI )X_ + 



Note that sin 0(a) = 

1 + cos 2 a 

obtain an exact expression for M(s max , a). 



Corollary C.2. We have 



and cos 9 (a) 



' -cos 9(a) -sin0(a^ 

1 

^0 — sin 0(a) cos 9(a) j 

( -cos 9(a) wx9(a> » 

1 

^0 sin 0(a) cos0(a)y 

3 cos 2 a — 1 



1 + cos 2 a 



By Lemma 



C.l 



we 



/ 



M(0,a) = M(s max ,a) 



1 

cos 40(a) sin 40 (a^ 
^0 -sin 40(a) cos 40(a) J 



(CI) 



where 9(a) is defined in Lemma C.l 



In other words, M(0, a) and M(s max , a) are rotations around x\— axis of angle 40(a). 
This fact is crucial for the derivation of suboptimal strategies presented in Sec. V A| It is 



worth noticing that formula (CI) is exact for any a smaller than 7r/4. If a is small enough, 
we have 40(a) = — 4y / 2a + o(a) which agrees the first order approximation used in the proof 



of Proposition IV. 1 



Corollary C.3. Starting from the north pole, the switching points of the extremals having 
their first switching at s max are located on the great circles passing through N and containing 
the x\— or %2—axis. 

Proof. Note that the switching points of these extremals are given by 



m 2 (a)M n (s mBX , a)m 1 (a)N, 
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where n is an integer, mi (a) denotes one of four exponentials in Lemma C.l, and 7712(0;) : 



e SmaxX++ , or e s ma X x + _ e s max x ++ ^ Qr e s m axX — e s max x + _ e s max x ++ _ Corollary Qj j g then proved by 



induction on n, using Lemma C.l and Corollary C.2 



□ 



The following two lemmas are valid for a small enough and /3 = it /A. 

Lemma C.4. Let v(s) be the second switching time as a function of the first one s. For a 
small enough, we have 



71 



(C2) 



where s 



7T 1 1 
max = - + + Y2° 4 + ^ 5 )' a7ld 



fx(s) 



- + cos s + sin s, 



25 1 . 1 ■ 2 

ms) := sins H — coss + coss sin s — cos s. 

J K 1 24 3 6 



Lemma C.5. We have 
( 



M(s,a) 



1 + f 3 (s)a 4 



h(s)a 2 + f 5 (s)a 4 



s cr 



-f 4 (s)a 2 -f 5 (s)a\ 1 - 16a 2 + f 7 (s)a 4 -A^a - f 8 (s] 



a 



s or 



\ 

M a (s,a) + 0(a 5 ), 



Ay/2a + f 8 (s) 



a 



16a 2 + f 9 (s)a A j 



+ O(o £ 



where 



Ms) 
Ms) 

Ms) 
Ms) 



16 sin s — 16 cos s sin s — 16 + 16 coss, Ms) := 4 — 4 coss — 4 sins 

2 s, 
2^2 



70 58 64 4 2 

— + — coss + — sins + 4 cos s, J6(s) := 8v 2(— 1 + cos s + sin s), 



112 

IT 

160 



16 coss sins, 



16 sin s — 16 cos s. 



-(—34 + 3 cos s + 3 sin s), 



Moreover, M _1 (s, a) = M~ 1 (s, a) + 0(a 5 ) = Mj(s, a) + 0(a E 
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